Energy absorption of a Bose gas in a periodically modulated optical lattice 
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We compute the energy absorbed by a one dimensional system of cold bosonic atoms in an optical 
lattice subjected to lattice amplitude modulation periodic with time. We perform the calculation for 
the superfluid and the Mott insulator created by a weak lattice, and the Mott insulator in a strong 
lattice potential. For the latter case we show results for 3D systems as well. Our calculations, based 
on bosonization techniques and strong coupling methods, go beyond standard Bogoliubov theory. 
We show that the energy absorption rate exhibits distinctive features of low dimensional systems 
and Luttinger liquid physics. We compare our results with experiments and find good agreement. 



Cold atoms provide a remarkable laboratory to study 
the physics of strongly correlated quantum systems. Cold 
atomic gases loaded in optical lattices 0, |2 and Fesh- 
bach resonances Q allow for an unprecedented control 
of many parameters of the system, including the interac- 
tions, both for bosons and fermions. However, the ability 
to probe the properties of such systems, and in partic- 
ular, to measure momentum and frequency dependent 
correlation functions still remains very limited. One of 
the most common experimental probe is time-of-flight 
(TOF) imaging, which under most conditions give ac- 
cess to the the momentum distribution M. More re- 
cently other spectroscopies such as Bragg , energy 
absorption rate (EAR) 0, 0, radio frequency and 
shot-noise have been demonstrated. Thus in order 
to properly characterize these systems it is of primary 
importance to have a good understanding as well as the- 
oretical predictions for such experimental probes. 

In a series of recent experiments, the Zurich group 0,0 
measured the EAR of Bose gases of ultracold ^''Rb atoms 
loaded in optical lattices. The system was heated by 
periodically modulating the lattice along one direction, 
and the energy absorption rate as a function of the mod- 
ulation frequency was subsequently estimated from the 
width of the atom distribution around zero momentum 
measured in TOF. Interestingly, of all the system dimen- 
sionalities reported in [3], the one-dimensional (ID) ones 
exhibit the broadest continuous spectrum away from the 
Mott regime, which has no simple explanation in stan- 
dard Bogoliubov theory I2.ri0ll. and has led some to con- 
sider other schemes 0, ll2Lll3l IT^ . 

In this work, we show that treating the ID interact- 
ing Bose gas in the framework of the Luttinger liquid 
(LL) 0, yvdi [l8| correctly takes into account short 
distance correlations that are neglected in the Bogoliubov 
theory, and leads naturally to a continuum of excitations, 
hence a continuous absorption spectrum. Furthermore, 
these correlations are responsible for the transition to the 
Mott insulating (MI) phase for sufficiently strong inter- 
actions. 

In contrast to the broad spectra observed for the su- 
perfluid (SF), in Refs.i and It] two much narrower peaks 
were observed for the Mott phase. Here we also study this 



regime and consider two physically distinct situations: 
First, a sufficiently strongly interacting Bose gas in ID, 
as in Ref. 0, becomes Mott insulating for a weak com- 
mensurate periodic potential 0, |^ . In this case, 
we present results for the EAR spectrum obtained using 
the form factors for the effective low-energy sine-Gordon 
model 22]. These results also apply to the MI phase of 
the ID Bose-Hubbard model near the SF-MI transition. 
We find that the spectrum exhibits a finite threshold, 
corresponding precisely to the Mott gap. Second, we 
consider a very deep optical lattice, so that hopping is 
strongly suppressed and the Mott gap is large. This sit- 
uation is well described by the Bose-Hubbard model, for 
which we have computed in \D, the shape of the lowest 
excitation peak in the absorption spectrum, which occurs 
at an excitation frequency lo « U /h, U being the on-site 
boson-boson repulsion. We also give the results for the 
peak width in higher dimensions. 

Consider a system of interacting bosons in a lattice cre- 
ated by an optical potential V{x,y,z) — VoxSin^{kx) -\- 
VoySm^{ky) -\- VozSm^{kz), where the wave vector is 
k = 27t/X, a being the laser wavelength. We first con- 
sider an optical lattice that is very deep in the y and z 
direction {e.g. Vqx <C Vj_ = Voy — Voz). Atoms accu- 
mulate in the minima of this potential, where they ex- 
perience strong transverse confinement, thus forming ID 
gas tubes with a weaker periodic potential along the axis. 
For large enough Vby and Vqz each ID system becomes 
isolated from each other |2^: this is the ID limit. 

To obtain analytical results, we neglect trapping and 
finite size effects 35]. If Vqx <C /i, where fj, is the chemical 
potential of the ID interacting gas, the system in the 
presence of the lattice is well described by the following 
sine-Gordon (sG) model d El US : 
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where 4>{x) and Tl{x) are canonically conjugate fields; (fi 
represents density fiuctuations, and 9 — n dx'Il{x') 
corresponds to phase fluctuations. Vs is the speed of 
sound and the coupling to the axial potential is go ~ 
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PoVox{Vox / fJ,)"''~^ , with no the number of bosons per po- 
tential welL K is a, dimensionless parameter determined 
by the strength of the boson-boson interactions and the 
Unear density po fvA |23| : the SF-MI transition occurs at 
the universal value Kc = 2 (in terms of 7 = mg/}T?pQ, 
the dimensionless interaction strength of the Lieb-Liniger 
model with m the atom mass and g proportional to the 
scattering length, it corresponds to 7 ~ 3.5). 

A weak time dependent modulation of the lattice 0, 
Vqx Vqx + 6VxC0s{ujt), leads to a perturbation that 
can be written in the above low-energy description as 
H'{t) = 5Vx cos{ujt)0, with O = fo J dx cos(20), where 
/o = dgo/dVox- Within linear response theory, the EAR 
per particle at frequency uj is given by 
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where Xc'(^) is the Fourier transform of the retarded 
correlation function -ih-^e{t){[O{t),O{0)]) , with 6 the 
step function. Note that, in the bosonization tech- 
nique 0, 01, the density p{x) — pQ — dx(t>{x)/'K + 
cos [20(a:) — 2'Kpox\ and hence the EAR probes the 

q ~ 27rpo = 2nofc part of the excitation spectrum. Stan- 
dard Bogoliubov theory has no low energy excitation near 
this momentum, hence no absorption, and one has to re- 
sort to non-linear response to account for the observed 
spectrum But in ID, the correct excitation spec- 

trum has in fact a continuum at q ~ 27rpo (25i. .26.] , which 
is taken fully into account by the bosonization method. 

In the ID SF (namely, LL) phase, where K > 2, the 
cosine term of ^ is irrelevant and the system is gapless. 
Xo{'^) for fiw ^ ^ can be obtained [23 by means of 
bosonization techniques 0,0], and we get: 
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where A is a non-universal prefactor depending on the 
microscopic details of the model (Bose-Hubbard or Licb- 
Linigcr). Eq. directly shows that the continuum of 
low energy excitations at q » 27rpo leads to this con- 
tinuous absorption curve [Sg. As Huj ^ p, bosoniza- 
tion ceases to be valid. However, standard sum rules for 
the density response function |28l | imply that the inte- 
grated absorption spectrum is finite, and thus the spec- 
tral weight must decrease as w ~ p/h (or to ~ J/h for the 
Bose-Hubbard model). One thus expects a rather broad 
spectrum in ID, as observed in the experiments 0. 

On the other hand, in dimensions higher than one, 
continuous absorption must start from a finite frequency, 
as required by Landau's criterion for superfluidity. We 
speculate that the relatively broad absorption spectrum 
observed for the SF in the "ID to 3D crossover" regime 
may be due to the possible existence of a broad roton-like 
feature. Since this particular lattice is highly anisotropic, 
a naive extrapolation of the dimensional crossover theory 
of Ref. [23] suggests that at high frequencies, the ID-like 
excitation continuum persists, while at low frequencies. 
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FIG. 1: Energy absorption rate (EAR) for ID interacting 
bosons in a weak optical lattice, for the MI phase (with small 
Mott gap) with different values of the parameter K (in brack- 
ets are the corresponding values of 7 [l3|)- 



the excitation is more 3D-like, with a broad roton-like 
minimum being the vestige of the ID continuum at q ^ 
27: pq, but this minimum must not be at zero frequency 
to satisfy Landau's criterion (see also Ref. .29|). 

Next we turn to the MI phase in ID, which occurs for 
sufficiently repulsive interactions, K < 2. Now, the co- 
sine term in leads to an excitation gap, with gapped 
solitons and anti-solitons in the 1 < K < 2 range of in- 
terest here 0, . Using the form factor approach, and 
keeping only the one soliton- one antisoliton contribution 
to the absorption, we find j2^ : 
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where i? is a non-universal prefactor, A — 2AIgVg is 
the Mott gap, and Ms - p{go/ papf/'^'^"'^'^ /vl {K < 2) 
is the soliton gap (mass). 6o{oj) — 2 arccosh(?kj/A) 
is the relative rapidity of the soliton and the antisoli- 
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factor of the operator cos 20 22], with ^ = K/{2 — K). 
The EAR is plotted in Fig. ^ for various showing the 
Mott gap to energy absorption. For ^iw > A the EAR 
increases monotonically in a way slower than in the SF. 

We now consider the case where the lattice potential 
is very deep {p <C V^x), for a c?-dimensional hypercubic 
lattice. The system is then described 'soj by the Bose- 
Hubbard model -ffsH — Hj + Hjj , where 
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creates a boson at lattice site R, Xa are lattice vectors 
joining the site R to its two nearest neighbors along direc- 
tion a — 1, ... ,d; Srip^ — — no. For Vba ^ Epi lat- 
tices, J„ = {%/^)ER{Voo./ERfl^eM-'^{Voo.lERf/^] 
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and U = AV2^{as/\)ER{Vo^VQyVQjElf'^ with 
as the scattering length, can be controlled by varying 
the laser intensity Vqq (measured in units of the recoil 
energy = h^k"^ /2m). In one dimension {i.e. rf = 1, or 
effectively when Ji > , for a = 2^. . . ,d), the SF-MI 
transition occurs at {U/J)c — 1-92 [sj, while in a d = 3 
square lattice, [U/J)c — 5.8, with J = Jq. Near the 
transition, on the MI side, the above description in terms 
of the sG model (UJ, still applies Thus, the EAR is 
also given by Q), and the Mott gap A is exponentially 
small. As U / Ji grows, eventually A w [/. 

Since Ja = Ja[Voa) and U — U{{Voa}) are functions of 
the optical potential strength Vba , the modulation along 
one direction Vqx = Voi Vvi x + 5Vx co&iujt) induces the 
following perturbation to © |33 |: 

H'{t) = SV^Fu Hbh +OcosLut, (6) 

where O = -^Y.r SJ^ ^^^^+3,^ , with SJi = Ji{Fj - 
Fu)5Vx, 5 Jo, = - jJFuSV, (a > 1), Fj = dlnJ^/dVo, 
and Fjj = d\nU /dVox . The first term in © is cx Hbh 
and does not contribute to the absorption. 
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FIG. 3: Width at the base of the first resonance peak in the 
MI phase. In the ID, 1D-3D crossover and 3D cases we take 
V± — 30Er, V± = 20Er and V± = Vxo respectively, and 
no — 1. The 3D Exp. case is the same than 3D but with 
no = 2, which is closer to the experimental values. The inset 
compares the predicted half widths with the experimental 
data from Ref.'^for ID. 
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FIG. 2: Spectral weight of the first resonance peak in the MI 
phase in ID. In the inset we show the predicted ID shape of 
the peak for U/ J = 36, and the comparison with (normalized 
in the vertical axis) experimental results from Ref. ^ 



The EAR can thus be computed by studying the lin- 
ear response to the last term of via a strong coupling 
expansion in Ja/U. We sketch here the derivation in the 
ID case, full details and the more involved higher dimen- 
sional case can be found in [23|. In the MI phase with 
uq bosons per site {uq is an integer), the ground state 
of Hij is |<I>o) = \ni — riQ, . . . , um = no), while the first 
excited state |$(i?,r)) with the same number of bosons 
has an extra boson ( "particle" ) at site R = 1, . . . ,M and 
one fewer boson ("hole") at R + r (r — I, . . . , M — 1) 
and costs an energy U. To obtain the linear response 
for hiu !v U, we only need the ground state and the 
states with one particle and one hole. To take into ac- 
count one needs to diagonalize the kinetic energy 
in a subspace where the particle and the hole hop, but 



no additional pair is created or destroyed. Thus, at 
arbitrary filling no, we employ the ansatz \^{Q,q)) = 

#Efl,re*'^'^e'^^'^^''sin(f ) |$(i?,r)) (Q = 2ttj/M and 
q = 2^1 /M for M lattice sites with periodic boundary 
conditions). The eigenvalues of the kinetic energy are 
e{Q,q) = [/ - Jp{Q)co8{q/2). Here p{Q) and e{Q) are 
the modulus and the argument of tiq + 6**^(1 -I- no)- Us- 
ing the spectral decomposition of the correlation function 
Xq{i-^) (to leading order in perturbation theory in 6Ji), 
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we get the EAR in ID to be: 
SJf / no + 1 
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There is thus a resonance at huj ^ U with a width at 
the base of the peak 2W = 2(2no -I- l)Ji. At lowest 
order in SJi the absorption is zero for fko < U — W or 
hu! > U + W. Note that the shape of the resonance peak 
is not symmetric around Hcv = U. 

This is shown in the inset of Fig. 13 where the ab- 
sorption maximum is at ^iw « U[l + (2no -I- 1)^J^/C/^]. 
The square root form in Eq. ((SJ is a non-perturbative 
result, i.e., an effective resumation of a certain class of 
diagrams [23. The peak spectral weight, /, can readily 
be computed: 
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and is shown in Fig. [3 In d > 1 there is no simple formula 
for the precise shape of the resonance peak. However, the 
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width of the peak still has the same form, provided that 
Ji is replaced by J. Moreover, it does not depend on (5Jq, 
and therefore it is unaffected by the a lattice modulation 
in one or more dimensions. In Fig. |3| we plot 21^ as a 
function of U / J, for the same setup as in Ref.|2i. The line 
shape shown in the inset of Fig. 3 is in good agreement 
with the experimental observations reported in Ref. |2, 0- 
Finally we compare our results with the experimental 
findings. In the strong coupling MI phase, our results are 
in good agreement with the experimental observations: 
(i) In ID, the line shape shown in the inset of Fig.|21is in 
good agreement with the experimental findings reported 
in Refs. and[l (ii) As shown in Fig.^l for U/J > 10, 
the width of the peak at U decreases with U/J in all 
dimensions, which is again in agreement with the exper- 
imental observation. Our values for the half width in ID 
are also in good agreement with experiments, as shown 
in the inset, except for a slight broadening that is due to 
the trapping, and also to the fact that the experiments 
are done in a multiple tubes setup, which changes the 
filling from one tube to another, (iii) In experiments, for 
fixed U/ J, the width of the resonance becomes smaller 
as one goes from the ID to 3D lattices. Our calculations 



capture this effect in the ID and 1D-3D crossover cases, 
although a larger value of the width at the base is pre- 
dicted for 3D [33. However we would like to emphasize 
that here we compute the width at the base whereas the 
half width is fitted in the experiments. Clearly, more 
experimental data in the MI regime are needed to test 
not only for the width of the resonance peak '3^, but 
also its spectral weight. In the weak coupling regime, as 
noted above, our results are also able to explain qual- 
itatively the broad absorption spectra measured in the 
experiments for ID systems, with a continuous absorp- 
tion curve for the SF, starting from uj = 0. Accurate 
measurements of the low-frequency part of the spectrum 
should allow, in principle, to determine the parameter 
which characterizes long-range correlations in the LL. 
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